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Abstract. We propose a new stochastic method to describe low- lying excited states of finite 
nuclei superposing multiple Slater determinants without assuming generator coordinates a 
priori. We examine accuracy of our method by using simple BKN interaction. 



1. Introduction 

Self-consistent mean-field models using effective interactions (Skyrme, Gogny, etc) have been 
applied to description of ground-state properties of nuclei with great success for a wide mass 
region^ |2]. In addition, the generator coordinate method (GCM) is often utilized to take into 
account configuration mixing. However, the GCM requires us to choose generator coordinates 
with physical intuition and the coordinate is restricted to one dimension in most cases. Thus, 
it is difficult to describe many kinds of low-lying excitations simultaneously. To overcome these 
difficulties, we attempt a new approach to select many Slater determinants. 

2. Method and results 

We briefly discuss how to select Slater determinants \§ s ) (s = 1,2, ■■■,N). We use the 
imaginary-time method which is a standard procedure for solving self-consistent Hartree-Fock 
equations. The energy expectation value of a Slater determinant is lowered as the imaginary 
time evolves and finally converges to the energy minimum state. If the configuration approaches 
to a local minimum or a nearly stable quasi-minimum (shoulder state), the state stays unchanged 
for a long period of time (FigQ). We would like to select Slater determinants corresponding to 
these local minima and soft modes near the ground state. For this purpose, we start imaginary- 
time calculations from many initial Slater determinants which are randomly generated, then 
pick up a state every 250 iterations of the imaginary-time step. Since the state stays near local 
minima and the ground state for a long time, this prescription leads to a large probability of 
selecting states we want. Repeating the same calculation with many initial states, we may select 
large number of Slater determinants without assuming generator coordinates. We project each 
state on parity and angular momentum [3], then perform the configuration mixing calculation 
by solving the Hill- Wheeler- type equation 
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Figure 1. An example of the path generated 
by the imaginary-time calculation. Energy as 
a function of the iteration and the density 
plots at every 500 iterations are shown. 



-100 

-105 - 

-110 - 

-115 - 

-120 - 

-125 - 

-130 - 

-135 : 

-140 - 




4 6 8 10 

Number of Basis States 



Figure 2. Energy of J w = + states in 16 
as a function of number of basis states. Solid 
and dotted lines indicate two independent 
calculations. 



It should be noted that, in order to avoid overcompleteness of selected set of basis, we exclude 
a state |3>) if this does not satisfy the condition, 

m^R n \^}\<0.7, (1) 

for any of selected states {I*!* 5 }}. R n (n = 1, 2, • • • , 24) indicate special rotations corresponding 
to permutation of the axis (x, y, z). P is the parity-inversion operator. 

Let us demonstrate that our method produces identical results for different initial sets of 
configurations generated by independent random numbers. In Fig|2J energy of J w = + states 
in 16 are shown as increasing number of Slater determinants. Solid and dotted lines are 
independent calculations starting from different initial states. We see that the calculated energies 
of lowest four levels are almost identical in these two calculations. This indicates a qualitative 
convergence of the result with respect to the number of basis states. Note that the obtained 
ground-state energy is lower than that of a single Slater determinant by about 2 MeV. 

3. Conclusion 

We proposed a new stochastic method for description of low-lying exited states of finite nuclei 
without assuming generator coordinates. Each state is composed of superposition of multiple 
Slater determinants to include the correlation beyond the mean-field. Basis states are generated 
by imaginary-time evolution with randomly-generated initial configurations. The method is 
suitable for selecting local minima and soft modes. We tested accuracy of our method by using 
simple BKN interaction. Calculation with more realistic effective interaction is a future subject. 
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